Quantum Heisenberg Antiferromagnets 
versus Nonlinear a Model Without the Large S Limit 
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In this letter, I develop a new topologically invariant coherent state path integral for spin systems, 
and apply it to the quantum Heisenberg model on a square lattice. As a result, the quantum 
nonlinear a model for arbitrary values of spin can be directly obtained. The effective coupling 

constant and spin wave velocity are modified by gs — 



Ta 



d + 257 and c. 



Ta is a natural temperature scale for the reliability of the theory, 
extended to other generalized coherent state path integrals. 
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The formulation can also be 



d + 257 , where 



Physics of quantum Heisenberg antiferromagnets 
(QHA) in low dimensional strongly correlated systems 
continuously attract attentions. This is largely due to 
the intense interests in understanding high Tc copper- 
oxide superconductivity which arises as a consequence 
of hole (or electron) dopings from the parent antiferro- 
magnetic (AF) Mott insulators 0. The AF Mott insula- 
tors concerned in high Tc basically correspond to a two- 
dimensional spin-i lattice QHA with nearest-neighbor 
exchange interaction, which is, quantal mechanically, still 
a difficult problem to solve. Meanwhile, there are also 
increasing interests in quasi-one-dimensional Heisenberg 
spin chains because of the stripe phase evidenced recently 
from the observed incommensurate magnetic excitations 
in cuprate superconductors l2l . In the large spin S limit, 
Haldane has shown that the lattice QHA can be described 
by an quantum nonlinear a model (QNLcrM) when the 
correlation length is vary large, i.e. at low temperatures 
p|. Since then, the QNLcrM has become a good candi- 
date to describe phenomenologically various experimental 
data of the two-dimensional QHA [0-pl . 

However, two crucial questions arise for some time: 
why the predictions of the long- wavelength, low energy 
physics of the QNLcrM obtained in large S limit can tru- 
ely be identical to experimental data of the lattice QHA 
in vary small values of 5'? and which is the actual tem- 
perature upper limit where the QNLcrM is reliable? In 
fact, it has been found that the correlation length £,{T) 
derived from the QNLcrM |^ is not always agreement 
with the experimental data of the QHA 0. Also, the 
temperature dependence of £,{T) does not occur the so- 
called quantum critical regime predicted by QNLcrM [gj . 
In this letter, by improving the formalism of spin coher- 
ent state path integral, I derive the QNLcrM form the 
lattice QHA for arbitrary values of spin. Meanwhile, a 
temperature scale Ta for the reliability of the theory is 
naturally realized. 

The general spin coherent state path integral is known 
for a long time 19] . Explicitly, the partition function can 
be expressed (in units ?i = fc^ = 1) by 



Z = Tre 
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where the term iS j^ drA <l = j^ dT{n\^\n) = 
iSuj{ft) is a topological Berry phase at the site i [ [10[ , 
and A is a U(l) monopole potential, the state |rj) is a 
spin coherent state which can be constructed by rotating 
the "north pole" state 15*, 5) to along the ft direction, 
\ft) — R{ft)\S, S). In other words, SI is a unit vector 
along which the spin operator with quantum number S 
is maximally aligned in \ft). For the lattice Heisenberg 
model, H = J X]<ii> Si • Sj {J > 0) which has the global 
SO (3) spin rotational symmetry, one can obtain 



Zh = / v[n,] 



J'^dT-^iS^. AiJ7i-JS'^^ .. fJi-nj ? 



(2) 



The whole exponent in Eq. (0) represents an effective 
action of the Heisenberg model in terms of the spin 
coherent state path integral. By minimizing H{n) = 
JS^ S<ii> ^i'^j , one can find the classical ground state 
(Neel state) which spontaneously breaks the SO (3) sym- 
metry. Then expanding the action around the ground 
state, one can easily derive the spin-wave theory for the 
QHA that describes the long wavelength spin modes Q . 
However, the spin wave theory is only applicable for 
long-range ordered phase where the SO (3) symmetry is 
spontaneously broken. Based on Mermin- Wagner's the- 
orem P4 I , no symmetry can be spontaneously broken for 
T > in one and two dimensional Heisenberg model. In 
other words, the spin wave theory is invalid in low dimen- 
sions. To derive an effective long wavelength action that 
keeps the full spin rotational symmetry, Haldane consid- 
ered the large S limit. In the large S limit, the path 
integrals of Eq. (p|) are dominated by the semiclassical 



equation: iSTi x i7 



dH[n] 

an 



By separating the semiclas- 



sical solution Cli into a slowly varying Neel order unit 
vector (— l)*n(a;i) plus a slowly varying magnetization 
density field perpendicular to n{xi) (Haldane mapping), 
then taking the continuous limit and integrating out the 



magnetic density field, Haldane shows that Eq. (g) can 
be reduced to the QNLcrM, 



Z}i ex I Vln] e 



ln]^-^Jd^+'.a,nO-n^ ^3^ 



which is defined in the d+l-dimensional space (x^, 
• • • , x'^^^)={x^ , • • • , x'', Cot), where d represents the d- 
dimensional lattice space of the Heisenberg model, A — 
a~^ is an intermediate momentum cutoff as the inverse 
of the lattice spacing a, go = Coa^^'^/p° — 2^/d/S the 
dimensionless coupling constant, Cq = 2ydJsa the spin 
wave velocity, and p° = JS^a?~'^ the spin stiffness. The 
imaginary time (temperature) variable r ranges from 
to /3 = 1/T. The exponent 8[n] in Eq. (H) is a topologi- 
cal factor associated with the Berry phase, which distin- 
guishes between integer and half-integer spins. 

In this letter, by improving the formulation (nl) of the 
spin coherent state path integral, I can indeed directly 
obtain the QNLcrM by integrating out the quantum fluc- 
tuation above the energy scale fc^TX. In this derivation, 
I neither assume the large S limit nor use Haldane's map- 
ping so that the result is applicable to arbitrary values 
of spin. Meanwhile, it also solves a long-standing prob- 
lem in the construction of generalized coherent state path 
integrals |14| . 

Simply speaking, the long-standing problem in the con- 
struction of coherent state path integrals arises from the 
fact that Eq. (|l|) is not well defined, so does Eq. (g). The 
main suspect comes from the assumption that |r2(T + 
5t)) — \^{t)) is order of 0{5t) used in the derivation 
of Eq. (|^). Although this assumption has been widely 
accepted in all the applications of generalized coherent 
state path integrals, it has never been justified [O. As 
Klauder first pointed out |l3], the semiclassical (or sta- 
tionary phase) approximation of Eq. (nj) [or (0) that Hal- 
dane used to derive Eq. (H)] is indeed incompatible with 
the required initial and final boundary conditions embed- 
ded in the coherent state path integral. 

Explicitly, coherent state path integrals are defined on 
a geometric space G/H (here is SU{2)/U{1) ^ S^) which 
has a phase space structure that the curvature du;(ft) 
obtained from the Berry phase serves as a symplectic 
structure for the corresponding semiclassical motion [g| . 
Then, for each pair {p, q) which obeys two first-order time 
differential equations, there exist two pair (i.e. four) 
boundary conditions (poi^o) a-nd {pf,qf) from the ini- 
tial and final fixed points, which is incompatible. In the 
literature, there are two possible ways to avoid this in- 
consistency. One approach was pointed out by Faddeev 
|l6| that one must specify the initial and final boundary 
conditions only by the independent complex variable Zi 
and z*r respectively due to the complex structure of quan- 
tum state M . Another approach is proposed by Klauder 
that one may add an additional square term of time- 
derivative, {tlY into the action of Eq. (nl) to match the 



additional initial and final boundary conditions in phase 
space p5| . 

Physically, the above inconsistency arises essentially 
from the ignoring quantum fluctuations in the coherent 
state action. In the path integral formulation, there al- 
ways exist simultaneously fast and slow varying paths 
that are associated with short and long range quantum 
fluctuations, respectively. The effective action for slow 
varying motions can be properly obtained by integrating 
over short range quantum fluctuations. However, in the 
derivation of Eq. (Q) , only these slowly varying motions 
are concerned and short range quantum fluctuations are 
simply dropped from the off-diagonal matrix elements. 
This causes the incompatibility between the semiclassi- 
cal approximation and the boundary conditions. 

To see more explicitly, we begin with the discrete form 
of the partition function derived exactly from the coher- 
ent state representation ^ 

N N 

Z= lim Y\dfi{Vl^)eyi^{^\n{Vl^\n^-^) 
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it}' (4) 



where 1^''^) = \^^} because of the periodicity of the 
trace), e = P/N is infinitesimal as A'' ^ oo. The assump- 
tion that irj*^) — jrj*''^^) is order of 0(e) only extracts 
slowly varying motions, where rapidly varying motions 
are neglected. As a result, the time continuous limit of 
(0) is just the conventional spin coherent state path in- 
tegral given by Eq. (|l|). 

To include the contribution of rapidly varying paths, 
one should expand the near-by coherent state overlap to 
the second order terms that either are exclusively slowly 
varying or include at least one rapidly vary term. The 
topologically invariant terms of such contributions can 
be uniquely expressed by 

lim Vln(rjfe|J7fc_i) = S* / drliA-fl 

fe=l •^" 

+ in-{'tix&n)- —5n-5n + ---\, (5) 

where the time derivative H and the variation (5J1 repre- 
sent slowly varying motions and short range fluctuations 
(i.e., rapidly varying motions), respectively. The param- 
eter TA is an intrinsic shortest timescale to distinguish be- 
tween slowly varying and rapidly varying motions, which 
I will discuss in details later. The second and third terms 
in Eq. (|5|) are usually neglected in the conventional treat- 
ment of path integrals. For the Hamiltonian term, since 
it is already proportional to e, I only keep the off-diagonal 
expansion up to the quadratic order of 50,: 
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where a, a' are indices of spin components. Substituting 
Eqs. (|,|) into Eq. (|), one has 

Z = lV[fl]V[Sn] cxp / drUsA ■ Cl - H[fl] 

dH[n]- 



isn X n 



an 



sn 






d^H[n] 



dfladfla' 



SnaS^a' 



}, (7) 



which describes both the slowly varying motion ft and 
the short range fluctuations 5fl. 

If one artificially drops the dynamical dependence of 
short range fluctuations [namely, the higher order terms 
originated from the off-diagonal matrix element of the 
Hamiltonian in Eq. (|6|)], and integrate out the remaining 
sn terms in Eq. ([7|). it turns out that: 

V[n] exp /" dThsA-tl-H[n]--STACt^\. 

(8) 

This is just Klauder's nonconvcntional coherent state 
path integral [15|. As one can see this term is origi- 
nated of integrating over short range quantum fluctua- 
tions. Note that the shortest timescale t\ — i^. It 
shows that the effective contribution of short range fluc- 
tuations is obtained by integrating over quantum dynam- 
ics above the temperature (energy) Ta. When T\ —> oo, 
this term vanishes. Namely, no short-range fluctuations 
are renormalized, as one expected |0]. Hence, the coeffi- 
cient of this topologically invariant metric term, i.e. Ta, 
serves as a factorization temperature scale between long 
and short range quantum fiuctuations. This temperature 
scale (or time scale in real-time dynamics) must be deter- 
mined from the dynamics of the original Hamiltonian. In 
other words, naively ignoring the dynamical dependence 
of short range fluctuations in (|7|) is physically not mean- 
ingful. Therefore, only in the absence of a Hamiltonian, 
Klauder's formulation is correct. 

On the other hand, it is interesting to see that if one 

requires vanish of the second term, iSfl xfl ^^ = 0, 

Eq. m) is just a variation expansion of (El) with respect 
to the semiclassical dynamics Haldane used M. How- 
ever, Eq. (M) here is derived by carefully treating the 
off-diagonal matrix elements of near-by coherent states 
in Eq. (H) in terms of the short range fluctuation Sfl and 
the slowly varying motions fl. There is no semiclassical 
expansion beginning with. Hence, it has no necessary 
requirement to take the semiclassical limit by letting the 
second term vanish. Instead of, one shall integrate out 
the short range fluctuations SO, to obtain a low energy ef- 
fective action that can describe the long wavelength spin 



modes. Since no semiclassical approximation is made, 
the resulting low energy effective action is valid for arbi- 
trary values of spin. 

Now, we can apply the new formulation of spin coher- 
ent state path integral (^ to the square lattice QHA. 
To specify the AF ordering, let the slowly varying Hi = 
{—iyn{xi), here the Neel order n{xi) is a unit vector 
|n(xi)| = 1. Then taking the space continuous hmit 
^i ^ a^ f ^''^ where a is the lattice spacing: 



H[n] = js^ Y^ n, ■ Oj 



<y> 



'dJS^N 



JS^ 



dH[n] 
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sn^o, 

6fla6flf3 



d'^xY.idknix) ■ dkn{x)l (9) 



fe=i 



2dJS^ 



(10) 

on^anr— a^ jd^^snix).snix). (ii) 

Substituting (|§-|ll]) into (0) and integrating out the short 
range fluctuation 5J1, I obtain, 

Zh oc fv[n] e'2-se[n] ^^p | _ ^ | " ^^^^^^ 
d'^x 






dn 



dr 



V,n|2]}. (12) 



This is the QNLcrM for low energy QHA with arbitrary 
values of spin. Note that two basic parameters, the di- 
mensionless coupling constant gs and the spin wave ve- 
locity Cs, in (n2h crucially depend on S and Ta: 



9s 



S 



Ta 
2SJ 



^2JSc 



Ta 
2SJ 



(13) 



While, the topological phase factor 27rS'0[n] = 
27r5^j(— l)'A(a;i) •n(xi) keeps the same as in Haldane's 
derivation. 

The difference between (H) and ( |l^ ) mainly comes from 
the contribution of the pure quantum fluctuation, i.e., 
the 1/ta term in (M) which is ignored in Haldane's map- 
ping for the large S limit ra but it plays an important as 
argued by Klauder jTsj. Also, this difference is indeed as- 
sociated with the temperature scale for the reliability of 
the QNLctM. Usually one thinks that there should be no 
intrinsic cutoff for the imaginary time variable r because 
quantum ffuctuations exist on all time scale in path inte- 
grals M. But a low energy effective theory constructed 
from path integral is defined by integrating over high 
energy dynamics, i.e., high energy quantum fluctuations 
above some energy scale Ta. Without such an intrinsic 
cutoff, namely, let 1/ta — > or Ta -^ oo, Eq. ([l^ ) is 
reduced to 

ZHO,Jv[n] e^^-^^H e^p{-^Jj'x |V.np]}. 

(14) 



where ps = JS'^a'^~'^ is the spin stiffness. Except for the 
topological phase, this is just the classical d-dimensional 
NLcrM rather than Haldane's d + 1-dimensional NLtrM 
P|. That is, without using this intrinsic short-time cutoff 
or naively ignore this cutoff, it should be hard to correctly 
extract the fundamental quantum effect of QHA. 

On the other hand, the lattice spacing a indicates 
the existence of an intrinsic momentum cutoff A in the 
d-dimensional momentum space: A = 2y/T:\r{d/2 + 
l)]^''^/a = L/a. Correspondingly, there must exist an 
intrinsic energy cutoff Ta = CsA/27r [|ll|. Combining 
with (fsl), it turns out that 



J 



SI? 

47r2 
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l&TT^d 
~1?~ 



(15) 



For d ~ 2 and S ~ 1/2, we have L ~ 2y^ and thus 
T\/ J ~ 0.97. This provides quantitatively a temperature 
upper limit for the reliability of the QNLcrM: 



< T/J < Ta/J ~ 1.0 



(16) 



Meanwhile, the spin wave velocity Cs can also be explic- 
itly determined from (O) and dlq). For La2Cu04 which 
is a typical d — 2 spin-1/2 QHA with a = 3.79A and 
J ~ IbQQK, I obtain (inset back the unit Ti) 



hcs — 2JSc 



Ta 
2SJ 



0.85 eV A 



(17) 



This is in excellent agreement with the experimental data 
ftcs = 0.85 ± 0.03 eV A |l|. 

Now one can see that the basic parameters in QNLcrM, 
i.e. the coupling constant gs and the spin wave velocity 
Cs, are unambiguously defined in terms of J,S,a and d 
in the spin system. Let yi — Xi/a{i — 1, • • • , rf), yi+i =^ 
Cst/q, Eq. (|l^) becomes 

T 2 - / r 

Zh oc jv[n] e'2-se[n] ^^p | _ _L /"* dy,+. 



29s 
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(18) 



which is a dimensionless QNLcrM, where L is only a ge- 
ometrical constant [see above Eq. (|l5|)]. The tempera- 
ture dependence of the coupling constant gs is given by 
the running coupling constant gs (T) which can be deter- 
mined by the standard renormalization group approach 
Mm . Further applications of our QNLcrM to thermody- 
namic properties of the QHA for T < Ta, i.e. the calcu- 
lations of internal energy, specific heat, correlation func- 
tions, staggered susceptibility, and correlation length, are 
straightforward for different values of spin and lattice 
dimension. The effect of topological phase that distin- 
guishes between integer and half-integer spins can also be 
studied. The results will be presented in separate pub- 
lications. Moreover, The construction of a low energy 



effective field theory from the extended coherent state 
path integrals developed in this letter is a general ap- 
proach, in which the use of the shortest timescale plays 
an important role in order to obtain a self-consistent ef- 
fective field theory. This approach can also be applied 
to other generalized coherent state path integrals |^ for 
the study of strongly correlated or strongly interacting 
systems. 
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